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Abstract For technical applications thermoelectric materials with a high figure of 
merit are desirable, and strongly correlated electron systems are very promising 
in this respect. Since effects of bandstructure and electronic correlations play an 
important role for getting large figure of merits, the combination of local density 
approximation and dynamical mean field theory is an ideal tool for the computa- 
tional materials design of new thermoelectrics as well as to help us understand the 
mechanisms leading to large figures of merits in certain materials. This conference 
proceedings provides for a brief introduction to the method and reviews recent re- 
sults for LiRh204. 



1 Introduction 

Against the background of climate change and the present energy crisis, the quest 
for alternative, green energy sources is more urgent than ever. In this regard, ther- 
moelectric materials which transform waste heat (gradients) into electrical power 
through the Seebeck effect H] 121 are particularly appealing. However, due to a low 
efficiency we have not yet witnessed a wider technological application almost 200 
years after Seebeck's discovery. Instead, thermoelectric al applications are restricted 
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to niche markets such as radioisotope power systems for satelhtes [i!3|. A possible 
first major application is the exhaust heat of cars and trucks, as automobile compa- 
nies presently test thermoelectrical generators in prototypes Such efforts could 
be put on another level if novel materials with a higher figure of merit ZT, where 
Z is the power factor and T the temperature, and hence a higher efficiency, were 
available. Most present technical applications use semiconductors such as Bi2Te3 
111 where recently power factors Z considerably larger than 1 could be achieved 
through phonon Q and bandstructure engineering 1 6 1 . 

Very promising are novel materials on the basis of strongly correlated electron 
systems (SCES) which are at the core of the present conference proceedings. 
This class of materials is very diverse, ranging from metals to Kondo insulators and 
semiconductors, from dto f electron systems, from relative simple crystal structures 
such as FeSb2 |8 1 to most complex metallic cage compounds. 

Having such a wide field and the additional possibilities to nano- and heterostruc- 
ture these systems, a better theoretical understanding and reliable tools to compute 
thermoelectric properties quantitatively are mandatory. Theoretical physicists from 
the SCES community have analyzed thermoelectric materials mainly on a model 
level, i.e., on the basis of the Falikov-Kimball, Hubbard and periodic Anderson 
model 121 Uni, often employing dynamical mean field theory (DMFT) ifTTI [Til [T3]| . 
These calculations showed, among others, the importance of correlation-induced 
enhancements of the effective mass generating a high, but narrow density of states 
-or spectral function to be precise- close to -but not at- the Fermi level. As a con- 
sequence, the thermoelectric figure of merits can be strongly enhanced. On the other 
hand, theoreticians from the density functional theory (DFT) L14j community have 
been emphasizing the importance of a particularly high density of states (DOS) 
llTSl HU and of the large group velovities for certain shapes of the bandstructure 

uni. 

Since both, correlations and bandstructure, can substantially contribute to en- 
hanced thermoelectrical figures of merit, we need to deal with both of them on an 
equal footing. Only if both aspects are optimized we can expect to design materials 
or artificial heterostructures with a really large figure of merits. Taking correlations 
and bandstructure into account is possible with the merger |[T8l[T9i of DFT in its lo- 
cal density approximation (LDA) li20J and DMFT, for which the name LDAh-DMFT 
was coined 121], see Il22l|23]l24l for reviews. While LDAh-DMFT has been applied 
already to many SCES materials, thermoelectrical properties have been calculated 
rarely in the past. Noteworthy exceptions are LaTiOa |25| and LiRh204 |26|. The 
main reasons for this is that a wider experimental interest in SCES thermoelectrics 
emerged rather recently and that the calculation of thermoelectric properties such 
as the Seebeck coefficient requires some additional postprocessing which is not yet 
standard in LDAh-DMFT calculations. 
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In the following, we will give a brief, elementary introduction to the LDA+DMFT 
approach in Sec. |2] This Section is divided into the three steps LDA (Sec. |2.1[ ), 
DMFT (Sec. |2.2| i, and the necessary postprocessing for calculating thermoelectri- 
cal response functions (Sec. |2.3| l. Sec. |3] presents exemplary results by hands of 



LiRh204 which are reproduced from Ref. |26|. Finally, Sec. 3.1 gives a summary 
and an outlook. 



2 LDA+DMFT method 

The aim of this section is to give the reader a brief, elementary introduction to the 
LDA+DMFT approach; for more details see the reviews ll22l l23l l24l . 

Starting point is the general ab-initio Hamiltonian for every material which, with- 
out relativistic corrections, reads in the Born-Oppenheimer approximation 



kinetic energy lattice potential Coulomb interaction 




It consists of three terms: L) The kinetic energy given by the Laplace operator 4,, 
Planck constant Ti, and mass nie for every electron /. 2.) The lattice potential given 
by the Coulomb interaction between (static) ions at position R/ with charge Z/e 
and electrons at position r, with charge — e. 3.) Finally, the Coulomb interaction 
between each pair of electrons / and j [note the factor 1/2 is needed since each 
pair is counted twice in Eq. ([T])]. Input for the LDA+DMFT calculation is usually 
the experimental crystal structure, i.e., the positions Ri as an adequate relaxation 
procedure to determine the /?,'s from theory still needs to be developed. 

While Hamiltonian ([T} is easy to write down, it is impossible to solve, even nu- 
merically, for more than ff{\0) electrons, since the movement of every electron 
is correlated with that of every other electron through the last term: the Coulomb 
interaction between the electrons. These electronic correlations play a particularly 
important role if electrons are confined in or d or /-electrons or in artificial nanos- 
tructures. For such systems the typical distance jr,- — ry| between two such electrons 
on the same lattice site (i.e., two electrons in the set of d- or /-orbitals around the 
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same ion) is small so that the Coulomb interaction and, hence, also the electronic 
correlations are strong. 



2.1 LDA step 

Since it is impossible to solve Hamiltonian ([T]l, we have to develop approximations, 
and arguably the most successful approximation so far are those developed within 
the DFT framework, particularly the LDA |20|. Strictly speaking, DFT only allows 
to calculate ground state energies and its derivatives but not bandstmctures and ther- 
moelectric transport functions. However, it turned out that the auxiliary Kohn-Sham 
Lagrange parameters gj. often also describe bandstructures very accurately, making 
bandstructure calculations one of the major applications of LDA. Interpreting the 
LDA Lagrange parameters as the physical (one-electron) excitation energies, i.e., 
the bandstructure, corresponds to replace Hamiltonian Q by the Kohn-Sham ||271 
LDA Hamiltonian 



LDA = 



h^Ai V- -e^ 1 r 1 



Inie Y 47reo Ir,- - R/ 1 J 47r£o |r,- - r 



(2) 

This Hamiltonian shows that the complicated electron-electron interaction causing 
the complicated electronic correlations has been replaced by two simpler terms: 
The Hartree term describing the Coulomb interaction of electron r, with the time- 
averaged mean density p(r) of all electrons and an additional term V}f'^ which 
aims at including the effects of correlations and interactions. 

However, the exact form of this term is unknown and certainly it is not local 
in r as approximated in the LDA. One can take the V^c of the jellium model ||28l 
which has a constant electron density and is only weakly correlated. Hence, it is 
not surprising that LDA bandstructure calculations fail for SCES |20|. For such 
materials, which are at the focus here, we need to take electronic correlations into 
account more profoundly. 

A possibility to do so is to take the LDA bandstructure of the less correlated 
orbitals but to supplement that of the more correlated d- or /-orbitals by explicitly 
taking into account the most, important local Coulomb interaction. This leads to the 
Hamiltonian 

=^ = L °^1/aCk,„<. + ^ L Uf,f nua ni„ra' -AeY^hi,,,^, (3) 

kIma ilama' ima 



Hlda 

where the first part is the same as the LDA Hamiltonian (j2]) but in second (instead 
of first) quantization and in k and orbital space (with I and m denoting two different 
orbitals) with creation and annihilation operators c^^^ and c^^^^, respectively. 
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Fig. 1 Illustration of the different elements of the Coulomb interaction matrix of Hamiltonian 
jsj. There is an inter-orbital Coulomb repulsion U', which is reduced by Hund's exchange / for 
a ferromagnetic spin alignment, and an intra-orbital interaction U. Orbital rotational symmetry 
relates these quantities as (7 = U' + 2J. 



The second term explicitly takes the local Coulomb interaction on the same ion 
site / into account. Typically only the Coulomb interactions for d (or /) / and m 
orbitals are considered here. These interactions are spin and orbital dependent be- 
cause of the the exchange matrix elements leading to Hund's rules, see Fig.[T]for an 
illustration. Let us note that in Hamiltonian (|3]l only the density-density terms are 
included since the inclusion of the spin-flip terms of Hund's exchange became only 
possible in quantum Monte Carlo (QMC) simulations lf29l with recent improve- 
ments ll30l[3T1[32ll33lCT. 

Finally the third Ae term subtracts those contributions of U already taken into 
account in the LDA to avoid a double counting. For a truly ab-initio calculation, 
U\ J, and Ae still need to be determined. To this end, screening has to be taken 
into account; and a possibility within the LDA framework is to employ constrained 
LDA, for details see ||24|. 



2.2 DMFTstep 

Having derived a multi-orbital many -body Hamiltonian ([3]l from the ab-initio Hamil- 
tonian ([TJ, we still need to solve it. A possible way to do so is to use Hartree-Fock, 
allowing for symmetry breaking with respect to the spin and orbital elements, i.e.. 

Homo' i I omo' i I omo' 

(4) 

where {hi„^(ji) is the average occupation of the orbital m on site / with spin a' . How- 
ever, in this LDAh-I/ ll35l approach electronic correlations are neglected through Eq. 
(j4]); and the only chance to reduce the Coulomb interaction energy is by a strong 
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Fig. 2 In DMFT, we approximate the material specific lattice Hamiltonian ([3J by a problem where 
the interaction is replaced by a self energy on all sites except for one. This DMFT single-site 
problem is equivalent to an Anderson impurity model which has to be solved self-consistently 
together with the k-integrated Dyson eq. ([6j. 



symmetry breaking. Hence, tendencies to magnetic or orbitally ordered phases are 
grossly overestimated, as is the tendency to open gaps. Even within these ordered 
phases many-body aspects such as spin-polarons are neglected as was shown in [|36l . 

A reliable approximation to include the local correlations induced by the local 
Coulomb interaction of Hamiltonian Q is possible with DMFT iXT] [12] [T3l. We 
cannot derive this approach in full detail here and refer the interested reader to lfT2l 
and ll24l . The basic idea is visuahzed in Fig.|2j We replace the local interaction on 
all sites but one by a self-energy E{co). This gives rise to an Anderson impurity 
model of a single interacting site in a medium %(fi)) given by the self energy and 
the interacting Green function G{(o): 

%i(o)-' ^G{co)-'+E{co) (5) 

This Anderson impurity model, defined by its non-interacting Green function % 
has to be solved self-consistently together with the k-integrated Dyson equation, 
where the LDA bandstructure e,^^'^(k) enters as a matrix in the orbital indices (Vbz 
denotes the volume of the Brillouin zone): 

Glico) = J^ [co+^-[e^^^-Ae]M-Ef,„ico)]-' . (6) 

From a diagrammatic point of view, DMFT corresponds to all (topologically dis- 
tinct) Feynman diagrams of which, however, only the local contribution for the self 
energy is taken into account. Hence, it is non-perturbative in the Coulomb inter- 
action but neglects non-local correlations between sites. Recent improvements of 
DMFT include such non-local correlations by taking a cluster of interacting sites 
instead of a single one in Fig. [2| l37l [38l 1391 or by extending the diagrammatic con- 
tributions in the dynamical vertex approximation (DFA) |40|, also see t4ni42ll43]| . 

What we still need to do is to solve the Anderson impurity model self-consistently, 
which for realistic multi-orbital calculations is typically done by quantum Monte 
Carlo simulations, different approaches are discussed in [241 . The standard result of 
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such a DMFT(QMC) calculation is the interacting local Green function G(/a)v) for 
imaginary (Matsubara) frequencies icOy or its Fourier transform, the imaginary time 
Green function G(t). But also various correlation functions and susceptibilities can 
be calculated. 



2.3 Calculation of thermoelectrical response functions 

Starting point for calculating transport properties is the Kubo formula. For thermo- 
electric materials the Seebeck coefficient 

\e\AQ 

is of particular importance. It is given by the constants Boltzmann kg, unit charge e 
and the ratio of two correlation functions, the current-current and the current-heat- 
current correlation function 

Ao = lim f Axe'^' {T,j{x)m) (8) 

iv^O IV Jo 

Ai = lim^ - dxe<'' {TrmJom (9) 



iv^O IV Jo 

in the static limit, i.e., frequency iv 0. Here, is Wick's time ordering operator; 

and jgix) are the current and heat-current operators respectively. Also relevant 
is the heat-current-heat-current correlation function 

which yields the electronic contribution to the thermal conductivity K similar as 
Aq does for the electrical conductivity a. Since the phononic contribution to the 
thermal conductivity is however typically much larger at room temperature and can 
be reduced by phonon engineering, we will not consider K in the following. Instead 
we will concentrate of the purely electronic contributions to the power factor 

i.e., on S and a. 

Diagrammatically, the correlation functions Eqs. ( 8|9|10 i correspond to Fig. [3] 



As indicated, the vertex F is usually not taken into account. In case of full orbital 
degeneracy (of the low energy orbitals), this holds exactly since one can show by a 
simple argument that vertex contributions are, for the local DMFT vertex, odd in k 
and hence their integrated contribution vanishes, see [44]. In the case of LiRh204 
where the three low energy orbitals ai^e very similar in energy and occupation, ne- 
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Fig. 3 Diagrammatic representation of tlie (heat-)current-(heat-)current correlation functions Eqs. 
(8|9|10| witli an incoming frequency (wiggled line) iv 0. For the current operator the wiggled 
line yields a factor (k), for the heat-current operator a factor ejl^'^'^v'* (k). The vertex F is typically 
neglected as indicated so that the calculation of the correlation function reduces to the simple 
bubble diagram of two (interacting) Green function, i.e., a factor G{k, (o) for each of the two 
double lines. 



glecting the vertex is still justified approximately. This allows us to calculate the 
bubble diagram for the (heat-)current-(heat-)current correlation functions A„, from 
the spectral function p{k,(o) = — l/n ImG(k, (o).ln the .x-direction, we obtain for 
diagram Fig. [3] 

A,n = 2%h r dco i^Tr[v^(k)p(k,a))v'^(k)p(k,a))] /(a))/(-a))(/3a))'". 

J-OO V ]^ 

Here v'k are in the general formalism the dipole matrix elements which we replaced 
approximately by the simpler group velocity obtained through the derivative of the 
dispersion relation: 

These are relatively easy to calculate from the LDA bandstructure. Note, here the 
quantities v^, e^'-'^, and p(k, co) are all matrices in the orbital indices. 

What we still need for calculating the DMFT (heat-)current-(heat-)current cor- 
relation functions, is the k-dependence of p{k,(o). In the DMFT self-consistency 
cycle, one calculates however only the local Green function 

G(«a;v) = ^£G(k,/a)v) (13) 

^ k 

at Matsubara frequencies /Oy. From the DMFT G{i(Ov) or its Fourier transform the 
imaginary time G(t), we can determine the optical and thermal conductivity as well 
as the Seebeck coefficient in some post processing steps: 

First, we need the self-energy for (real) frequencies. The standard procedure Ii46l 
to this end is first to analytically continue the Green function to real frequencies. 
This is done by the maximum entropy method f45l yielding ImG(a)) at real fre- 
quencies (0 from G(t). From this, in a second step, the full G{(o) is constructed 
by Kramers-Kronig transformation. Third, that self energy is determined which, 
if plugged into the k-integrated Dyson eq. (|6]), gives the Green function which is 
closest to the QMC-determined G{(0). Finally from the self energy r(ft)) for real 
frequencies, we can determine 
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or its imaginary part p(k, co). 

For the LiRh204 calculations presented below it turned out that this standard 
approach does not work so well because the high spectral weight close to the Fermi 
level makes the analytical calculation very sensitive to the statistical QMC error On 
the other hand, we only need E{co) at small frequencies, of the order k^T, if we 
are interested in thermodynamical responses to at static fields (as S and a). Hence 
in ll26 l. we did a Pade fit |47 1 to E{iC0v) which works rather well for not too large 
frequencies. Comparing it with a polynomial fit allowed us to estimate the error in 
E{(o), see Fig.|6]in Sec.|3]below. 

Let us note the connection to the Boltzmann approach. This is obtained for non- 



interacting electrons and a constant-T approximation, i.e., calculation Eq. ( 12 1 for a 



self energy E{co) = — // T. This reduces Eq. (12 1 to 



k 



dm 

de 



gLDA(k) 
kBT 



(14) 



Note that, in contrast to the thermal and electrical conductivity, T cancels in the 
Seebeck coefficient since we divide A i by Aq. Hence, the exact value of the difficult 
to determine relaxation time is not relevant, as long as it is constant. 

For a better understanding of the microscopic origin of a large thermopower, at 
least as far as trhwe bandstmcture effects are concerned, we can approximate the 



Boltzmann Eq. (14i by summing only the states in a window ±^3?" around the 



Fermi energy (indicated by the tilde below): 

^0 ~ t£^v1 + v'I ; Ai « t£^v1 - vl (15) 

Here and v| are the typical (averaged) velocities above and below the Fermi level, 
respectively. For the current-current correlation function Aq these two contributes 
have to be added, whereas they have to be subtracted for the heat-current-current 
correlation function Ai. The reason for the latter is that a quasiparticle above the 
Fermi level carries a positive energy contribution relative to the Fermi energy, while 
we have a negative energy-contribution for quasi-hole excitations below the Fermi 
level. For getting an (absolutely) large Seebeck coefficient we need a large A 1 rela- 
tive to Aq. Since Ai is the difference of the same (positive) contributions which are 
added Aq, this requires the minuend to be much smaller than the subtrahend in Eq. 
( [T5] l or vice versa. This is possibly if either (i) there are many more states below 
the Fermi level than above (or vice versa for a large negative S) or (ii) the group 
velocity above the Fermi level is much larger than v| (or vice versa). Optimal 
would indeed be a combination of both. The route (i) can be heavily effected by 
electronic correlations, e.g., if we have a sharp Kondo peak directly above or below 
the Fermi level, but also bandstructure effects play a role. In contrast for mechanism 
(ii) the LDA group velecities (or dipole matrix elements) enter so that electronic 
correlations are not direclty relevant. 
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Fig. 4 Crystal structure of LiRh204, made up from Li04 tetrahedra and the RhOj octahedra [re- 
produced from L48J ]. 

3 An example: LiRh204 

Let us here briefly review the calculation of the Seebeck coefficient for LiRh204 of 
Ref. ll26ll . This mixed-valent spinel, see Fig. |4] was most recently synthesized by 
Okamoto et al. BSl . It shows two structural phase transitions: cubic-to-tetragonal 
transition at 230K and tetragonal-to-orthorhombic transition at 170K. For the high- 
temperature cubic phase Okamoto et al. reported a thermopower as large as 80/i V/K 
at 800K, which for a metallic system, is quite exceptional. Together with Na^Co02 
|[49|, it shows that transition-metal-oxides are promising candidates for thermoelec- 
tric application since, even in the metallic phase, large power factors {S^a) are 
possible. Coneming NaxCo02, these experimental findings led to some "heated 
discussion" on the origin of the large Seebeck coefficient on the theoretical side 
ll50l [TTl ISTI . This makes a reliable ab-initio calculation which can put the theoreti- 
cal ideas on a more solid fundament mandatory. 

Starting from the experimental crystal structure Fig.|4] the first LDAh-DMFT step 
is the calculation of the LDA bandstructure. Our results, using linearized muffin tin 
orbnitals (LMTOs) |52 |, are shown in Fig. |5] (left panel; dashed line). We further 
simplified the LDA bandstructure by a Wannier projection|53| of the LMTO wave 
funtions onto the subspace of Bloch waves, which were in turn Fourier transformed 
to Wannier funtions, see ll54l for details. Here, we even model the LDA bandstruc- 
ture by a two-band model (solid line). 

The next step is a self-consistent DMFT calculation. To this end, quantum Monte- 
Carlo simulations were used as an impurity solver ll29l . The Coulomb interaction 
parameters were estimated as {U,U\J) = (3.1, 1.7,0.7)eV from [55| and tempera- 
tures /3 = 1 /ksT —30, 34, 40 eV^' were considered. From the imaginary QMC self 
energy we obtained the self energy on the real axis (Fig.|6]) through a Pade and poly- 
nomial (Taylor) fit. As one can see from a comparison of the two fits there is some 
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[states/Ry./atom] 



Fig. 5 Left panel: Band dispersion of the effective 3-orbital Hamiltonian (solid line) and total 
LMTO band structure (dashed line) of LiRh204. Right panel: partial aig and density of states 
for the model. LDA [reproduced from 1261 1. 

uncertainty, but absolute differences are small, i.e., of (9(0.01) eV. Nonetheless, we 
proceeded with both self energies to have an estimate of the error 

From the self energy. Fig. |6] we can estimate the quasiparticle weight Z = 
(1 — dReE /d(o)^^ and the effective mass enhancement m* /m = 1 /Z. This effec- 
tive mass enhancement is actually not very strong, i.e., « 40% for the band and 
« 30% for the aig band. This indicates that electronic correlations are only inter- 
mediately strong for this compound, even though it is a transition metal oxide. The 
reason for this is the mixed- valent nature of LiRh204 which puts the orbital occupa- 
tion far away from a (more strongly correlated) integer filling. A second noteworthy 
aspect, we can extract from the self energy is the strong frequency dependence and 
asymmetry of the imaginary part of the self energy. This poses the question whether 
a constant-ImE, i.e., a constant relaxation time T approach as in the much less in- 
volved Boltzmann approach, works. 

From the (two) self energy of Fig. |6] we calculated the Seebeck coefficient using 
the formulas of Sec. 12.31 As one can see there are some differences in the Seebeck 
coefficient for the Pade (x symbol) and polynomial fit (* symbol), giving us an 
estimate of the accuracy of our calculation. Both are in good agreement with the 
experimental values |48|. 

Besides the LDAh-DMFT study, we also performed calculations (i) putting a 
constant-T self energy into the equations of Sec. |2.3| (+ symbol) and (ii) using di- 
rectly the Boltzmann equation (solid line). As one can see in Fig. [6]both agree, as 
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Fig. 6 DMFT(QMC) self energy calculated by the Fade approximation (left) and a polynomial fit 
(right) [reproduced from Ii26i 1. 



one can expect from theoretical considerations; but it is a good test in an actual 
implementation as two completely different programs based on different equations 
were employed. The Boltzmann equation yields a slightly too large Seebeck coeffi- 
cient S, albeit it still agrees surprisingly well with experiment. The reason for this is 
that LiRh204 is not strongly correlated. Besides, the two and aig bands are not 
strongly shifted with respect to each other by electronic correlations and have a not 
too different self energy. Hence, we are not too far from a situation were the self en- 
ergy is orbital independent. In this case, the DMFT spectral function is just a more 
narrow (quasiparticle renormaUzed) version of the LDA DOS with the same height 
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Fig. 7 Thermopower calculated by the Boltzmann equation approach and the constant-T method 
as well as by LDA+DMFT, using both the Pade approximation and a polynomial fit for the self 
energy [reproduced from 1261 1. 



at the Fermi level. Because of this, relatively weak electronic correlations do not 
strongly affect the Seebeck coefficient. We hence attribute the differences between 
Boltzmann approach and LDA+DMFT to the non-constant and strongly asymmet- 
ric ImL. This means that, in contrast to the constant-T approximation, the actual 
life time of quasi-holes is longer than that for quasi-particles. Let us emphasize that 
electronic correlations play a much more prominent role in other transition metal 
oxides, so that for these the Boltzmann approach will fail. 

Being confident, that the Boltzmann approach roughly describes the Seebeck co- 
efficient of LiRh204, we analyze Eqs. ([TSji. To this end, we plot in Fig.|8]the group 
velocity along the indicated paths thought he Brillouin zone, within the energy win- 
dow of \e~Ef\ < 3kBT at T ~ 300K. Fig.[8](upper panel) shows that v| is consid- 
erably larger than v'| in large parts of the Brillouin zone, particularly around the K 
and W point. The reason for this difference is a particular shape of the bandstructure 
very similar to the ideas proposed in ifTTl for Na;cCo02. This pudding -mold type of 
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Fig. 8 Group velocity squared (v|) along different directions of the first Brillouin zone for Rh+-''' 
(LiRli204; upper panel) and Rh+-' '"* (electron-doped LiRh204; lower panel), k point above the 
Fermi energy Ef are shown in yellow, those below Ef in red (reproduced from \2(>i ). 



shape is sketched in the inset of Fig. [8] for the full bandstructure see Fig.|5] In con- 
trast to the one band situation in NavCo02 ifTTl . we have however a double pudding 
mold. For the lower band, the pudding-mold shape leads to a very flat bandstructure 
above the Fermi level, flatter than a simple maximum because of additional turning 
points and minima. Consequently the group velocity above the Fermi level is very 
small and the Seebeck coefficient largely positive. 

Having understood the bandstructure origin of the large thermopower in LiRh204, 
we are now in a position to identify routes to even further increase the thermopower. 
As the sketch in Fig.|8]suggests the upper pudding mold band does not strongly con- 
tribute to the Seebeck coefficient or thermopower. Here, the situation is more like in 
a standard metal with group velocities being large above and below the Fermi level 
(this is the region between F and X and between F and L point in the main panel). 
Therefore, positive and negative contributions to the Seebeck coefficient roughly 
cancel for this upper pudding mold band. 

We can improve the situation however by electron doping which shifts the Fermi 
level to higher energy. Then the situation becomes very much the same as for the 
lower pudding-mold band before and, at the same time, the lower pudding-mold 
band is still contributing with the same sign because there are states below the Fermi 
level but no states above. As one can see in the lower main panel of Fig. [8] doping 
by 0.42 electrons, i.e., for a valence Rh^^ *^^, indeed leads to a situation where only 
the squared group velocity below the Fermi energy is large. 

We further studied this idea by calculating the thermopower and the power fac- 
tor for various Rh valences, using the Boltzmann equation approach. Note that we 
assumed the electron doping not to affect the LDA bandstructure expect for a shift 
of the Fermi level. We also neglected the energy and filling dependence of T, which 
should be present and affect p and, hence, the power factor (albeit not S). Depending 
on how the electron-doping is realized T might change because of disorder effects. 
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Rh valence 

°3 3.1 3.2 3.3 3.4 3.5 
Rh valence 

Fig. 9 (Color online) Power factor (normalized by its value at Rh valence=+3.5 ) and thermopower 
(inset) as a function of the valence of Rh, calculated by the Boltzmann equation. 

Nonetheless, we expect the tendencies to hold also for the power factor in experi- 
ments electron-doping LiRh204. 

As one can see in Fig. [9] (inset) the Seebeck coefficient strongly increases with 
electron doping, i.e., with reducing the Rh valence towards 3. However if the Rh 
valence is 3, the band is completely occupied so that while the Seebeck coefficient 
is large, the conductivity a = l/p becomes small. Hence for the power factor S^O 
(Fig.|9]main panel), there is a trade-off between a larger Seebeck factor and a large 
resistivity if the valence goes towards 3. This trade-off leads to a maximum for the 
valence Rh+^ ''^, afore shown in Fig.^ 



3.1 Summary and outlook 

We presented a brief introduction to the LDAh-DMFT approach for the realistic cal- 
culation of thermoelectric properties, including bandstructure and electronic corre- 
lation effects. We have shown that the LDAh-DMFT results for LiRh204 well agree 
with experiment. Furthermore, we identified the origin of the large thermopower in 
this material to be a particular shape of the bandstructure of the form of a (double) 
pudding mold. Even larger thermopowers can be obtained if the material is electron- 
doped, according to our prediction. For the particular material LiRh204 the micro- 
scopic mechanism for the large thermopower is foremost the bandstructure since 
electronic correlations are not very strong (the effective mass enhancement is only 
40% and even less for the flig band). This shows the strength of LDAh-DMFT to 
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unbiasedly identify bandstructure effects as the origin of large thermopowers where 
this is appropriate and electronic correlations where these prevail. 

For getting the optimal thermoelectric material, hetero- or nanostructure we 
likely need both ingredients. First, a good bandstructure such as the pudding-mold 
form discussed in the present paper which due to dramatically different group ve- 
locities above and below the Fermi energy yields an extraordinarily large Seebeck 
coefficient. And second, correlation effects which result in asymmetrical, sharply 
peaked renormalized spectra in the vicinity of the Fermi level which enhance the 
Seebeck coefficient as well. With LDA+DMFT, we have an ideal tool to scan and 
design a wide range of potential SCES materials on a computer, providing experi- 
mental physicists and chemists with valuable hints on how to improve the thermo- 
electric figure of merit. In the exemplary case of LiRh204 this would be through 
electron-doping the material. 
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